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Abstract. Cosmic crystallography is based on the principle that peaks in the pair separation histogram (PSH) of objects in 
a catalogue should be induced by the high number of topologically lensed pairs that are separated by Clifford translations, in 
excess to "random" pairs of objects. Here we present modifications of this method that successively improve the signal-to-noise 
ratio by removing a large part of the noise and then false signals induced by selection effects. Given the transient nature of the 
most readily available tracer objects, active galactic nuclei (AGNs), the former is possible because a natural filter for removing 
many of the noise pairs is available: when counting pairs of objects in order to create PSHs, only those with nearly identical 
redshifts need to be counted. This redshift filter (a maximum value of Az/z = 0.005) was applied to a compilation of AGN 
catalogues. Further noise was removed by applying a second filter, a maximum angle A8 = 0.075 rad, and a minimum number 
of pairs «p a i rs = 3 to find each "bunch of pairs" (BoP) where the vectors (in Euclidean comoving space) defined by pairs are 
required to be nearly equal, whereas in the PSH only the lengths must be nearly equal. These filters reveal significant signals, 
which, however, are due to selection effects. A third filter, a minimum length ^ se ] ec = 150 h Mpc between the (parallel) 
vectors in a BoP, is found to effectively remove these selection effect pairs. After application of these successive filters, no 
significant topological signal was found. 
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1. Introduction 

General relativity (GR) relates the local geometry and con- 
tent of space-time in the limit as one approaches any given 
point. Together with assumptions about approximate local ho- 
mogeneity and local isotropy 1 , the equations of GR relate the 
metric properties of space (curvature) and the average material 
contents and dynamics (density, expansion rate, cosmological 
constant) of the Universe. 

GR is, however, only a local theory that does not constrain 
global geometry: the curvature of space does not tell us whether 
the Universe is finite or infinite in spatial volume (or mass). 
Schwarzschild (1900, 1998) pointed out a century ago that a 
flat, homogeneous universe can be finite if the topology of 
space is not trivial, i.e. if space is multiply connected. De Sitter 
(1917) and Friedmann (1923) also noted this, while Friedmann 
(1924) pointed out that even a universe of negative curvature 
could be homogeneous but finite. In contrast to local geome- 
try, which is constrained by GR, there is neither consensus on 
a physical theory that would predict what the topology of the 
Universe should be nor any theory that would relate that topol- 
ogy to dynamical properties of the matter distribution, even 



1 For comments on "local" versus "global" homogeneity and 
isotropy, see Roukema (2002a). 



though work on some elements of what might be needed has 
begun, e.g. Dowker & Surya (1998); Hawking (1984a, 1984b); 
Zeldovich & Grishchuk (1984). 

Nearly all methods for detecting the topology of the 
Universe rely on the principle of multiple topological imaging 
of either discrete objects or of the thermally radiating plasma 
seen in the cosmic microwave background radiation. This phe- 
nomenon is called 'topological lensing'. Detection of topolog- 
ical lensing would allow us to directly study the evolution of 
individual objects over a cosmological time-scale, as the same 
object would be visible at different cosmic epochs. 

Although the amount and the quality of extragalactic ob- 
servations available is growing constantly, and many differ- 
ent observational strategies are possible, detection of topolog- 
ical lensing has so far proved to be very difficult. Reviews of 
different strategies include, e.g. Uzan et al. (1999a); Luminet 
& Roukema (1999); Roukema (2002b); Reboucas & Gomero 
(2004). 

One of the methods, "cosmic crystallography" (Lachieze- 
Rey & Luminet 1995; Lehoucq et al. 1996), was explained and 
applied to some catalogues of observational data, but no signif- 
icant detection has so far been obtained. 

Is further application of this method still worth trying? 
Arguments in favour include: 
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(i) The microwave background observations of DASI, 
BOOMERanG, MAXIMA, DMR and CBI lead to Q t ot - 1 
(Sievers et al. 2003). This estimate has been strongly un- 
derpinned by the WMAP-based result: £l M = 1.02 + 0.02 
(Spergel et al. 2003). The fact that the Universe inside of 
the observable horizon appears to be nearly flat makes cos- 
mic crystallography a viable method. 

(ii) Active galactic nuclei (AGNs) have short lifetimes as ac- 
tive (easily visible) objects, leading to (approximately) 
removable noise in the pair histograms. In this work, a 
method of removing these pairs and improving the signal- 
to-noise ratio is introduced. 

The main argument against applying cosmic crystallogra- 
phy is that recent attempts at constraints using the WMAP mi- 
crowave background data have so far failed to detect matched 
circles for torus-like models (e.g. Cornish et al. 2004). Also, 
analyses of the WMAP data are accumulating that suggest 
that the global topology of the Universe may be that of the 
Poincare dodecahedral space (Luminet et al. 2003; Roukema 
et al. 2004). Clearly, at any time well past the quantum epoch, 
t » 10~ 43 s, the comoving spatial section of the Universe can- 
not simultaneously be both a Poincare dodecahedral space and 
a toroidal 3-manifold. 

These reasons are not, however, sufficient for choosing not 
to develop and apply several tests for such a fundamental prop- 
erty of the Universe by independent methods, using indepen- 
dent sets of objects. The cosmic microwave background (CMB) 
remains difficult to interpret and filter for foreground contam- 
ination and for the integrated Sachs-Wolfe Effect, and even 
the internal linear combination WMAP map (ILC) is found 
to include significant foreground components at low / val- 
ues according to some authors (Naselsky et al. 2003, 2004). 
Moreover, even if individual pixels in CMB data, rather than 
statistical properties, were perfectly resolved into their respec- 
tive cosmological components, it would still be useful to have 
physically independent checks on CMB results. 

In this article, the method of reducing the noise in the cos- 
mic crystallography method — in the case of a catalogue of 
AGNs — is presented in Sect. [2] This involves successive fil- 
ters, parametrised by: 

(i) Az — a maximum redshift difference criterion to remove 
many of the noise pairs in a pair catalogue 

(ii) AO, «p a j rs — to remove pairs which are least likely to form 
"spikes" of topological origin 

(iii) i se iec — to remove selection effect spikes which imitate 
topological spikes 

Results are presented in Sect. 13.21 and are summarised in 
Sect. @] Notation adopted here for local cosmological param- 
eters includes Q m (present matter density), Q.\ (cosmolog- 
ical constant and/or quintessence constant), w (quintessence 
parameter), and the Hubble constant parametrised as h = 
Z/o/lOOkms- 1 Mpc- 1 . 



2. The method 

2.1. The cosmic crystallography method and its 
variations 

The cosmic crystallography method was originally proposed by 
Lachieze-Rey & Luminet (1995) and presented in more detail 
in Lehoucq et al. (1996). In brief, the idea was the following. If 
the Universe (more specifically: the comoving spatial section of 
the Universe) is multiply connected, then each image of a given 
object is linked to each of the other images of the object by the 
holonomies of the space. Without knowing the curvature or the 
topology we do not know what these holonomies are, but we 
know that they are isometries. Therefore, any pair of images of 
a single object can be related by an isometry. At the separation 
distance corresponding to such an isometry, the probability of 
finding two objects separated by that distance should be higher 
than in the absence of the isometry. In other words, the spatial 
two-point auto-correlation function of a catalogue of objects 
should show excess correlations at certain pair separations. 

Lehoucq et al. (1996) presented this as a non-normalised 
auto-correlation function, simply plotting counts of pairs of 
objects and calling the function a "pair separation histogram" 
(PSH). They carried out several simulations starting from an 
artificial catalogue that consisted of a moderate number of ob- 
jects randomly distributed in a fundamental polyhedron with 
equal edge lengths. In such an idealised case, the excess cor- 
relations that gather pairs linked by holonomies stand out dra- 
matically above the background distribution of other pairs, and 
so are termed "spikes". Lehoucq et al. (1996) pointed out that 
when applied to real data (they used a catalogue of galaxy clus- 
ters), many factors start to play a destructive role, i.e. by smear- 
ing out such sharp spikes, which makes signatures of the pos- 
sible existence of ghost images difficult to detect. Uzan et al. 
(1999a) went on to apply the PSH test to a quasar catalogue 
and found no topological signal. 

An alternative proposal for collecting "generator pairs" (or 
"type II pairs"), as in cosmic crystallography, was the the 
search for local isometries, i.e. to group together "local pairs" 
(or "type I pairs"), or more generally, local n-tuplets (Roukema 
1996). In Uzan et al. (1999b), the special case of n-tuplets with 
n — 2, i.e. pairs, was discussed and extended. 

Uzan et al. (1999b) extended the 2-tuplet case to provide 
a single statistic, equivalent to the auto-correlation function of 
the spatial two-point auto-correlation function, and named this 
the "collecting-correlated-pair" (CCP) statistic, which should 
have a high value in the presence of multiple topological 
imaging in a population composed of good standard candles, 
whether space is hyperbolic, flat, or spherical. Since there 
should only be a significant signal for the correct values of 
the curvature parameters Q m (present matter density), £2a 
(cosmological constant and/or quintessence constant), and w 
(quintessence parameter), a signal could potentially be detected 
as a relative signal without having to calculate its absolute sig- 
nificance. 

Application of the CCP method to a quasar catalogue again 
gave no topological signature. The authors were limited by 
computer time in spanning the full parameter space with high 
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resolution in parameter values, but the principle is certainly 
valid. 

Another way to apply a similar test is for the case where: (1) 
the curvature parameters' values are known, (2) the Universe is 
negatively curved, and (3) a particular hypothesis is made re- 
garding the 3-manifoldof space. This is known as the "pullback 
method" (Fagundes & Gausmann 1999). 

While the PSH test is designed to detect topological lens- 
ing, it has been pointed out by Gomero et al. (2002) that (1) 
it is likely to miss some part of the real signal of topologically 
lensed pairs, and (2) it would be insufficient to determine which 
3-manifold we live in: 

(1) The spikes of topological origin in single PSH's are only 
due to translations of the covering group, whereas correla- 
tions due to the other (non-translational) isometries mani- 
fest as small deformations of the PSH of the universal cov- 
ering space. Hence, some part of the real signal would be 
missed (if a rotation in the covering space is needed when 
matching faces). 

(2) Distinct Euclidean manifolds that admit the same trans- 
lations of their covering space would present the same 
spike spectrum, despite being different; therefore, the set of 
topological spikes in the PSHs is not sufficient for distin- 
guishing these compact flat manifolds, making it clear that 
even if the comoving spatial section of Universe is globally 
flat and observably multiply connected, the spike spectrum 
may not be enough for determining its global shape. 

Another way of stating (1) is that among the locally ho- 
mogeneous, flat 3-manifolds, the PSH method works for the 
3-torus, because the distance from a source to its nearest im- 
age does not depend on its location, but for all other locally 
homogeneous, flat 3-manifolds, the PSH method detects only a 
subgroup of the full holonomy group, i.e. the subgroup of pure 
translations (Weeks et al. 2003). 

However, neither of these limitations suggest that the 
method should not be used; any signal that appeared to be sig- 
nificant would imply a family of hypotheses to be further tested 
(depending on telescope time allocation). 

2.2. Modifications to the cosmic crystallography 
method 

Although the active phase of a typical galaxy is believed to be 
a recurrent phenomenon rather than happening just once in its 
lifetime (see e.g. Hatziminaoglou et al. 2001), the fraction of 
time during which a galaxy is in an active phase is likely to 
be small. For quasars, for which the comoving space number 
density at z ~ 1 — 2 is about a thousand times less than that of 
galaxies, the fraction of time during which they are active is 
likely to be similar, about 10~ 3 . This makes the chance of see- 
ing both members of a pair of topological images small, though 
not zero. Neglecting the "ephemeric" nature of quasars is the 
factor that may have led Uzan et al. (1999b, 1999a) to have 
missed a true signal in the catalogue they analysed. 



2.2.1. Maximum redshift difference: Az 

The short lifetime of AGNs implies that any pair of topological 
images of a single AGN is most likely to be seen during a single 
period of activity of the AGN, i.e. during a short time or redshift 
interval, Az. From the light travel time argument and assuming 
that the length of an activity period is of the order of 10 7 years, 
we adopted Az/z = 0.005 as a reasonable amount for a wide 
range of redshifts. 

By removing all pairs of objects at zi,Zj for which 
\Zi - Zj\/Zi > Az/z, the number of potential topological image 
pairs in the catalogue will be unchanged, while the total num- 
ber of pairs will decrease by a large factor. Hence, the signal- 
to-noise ratio should increase noticeably. 

2.2.2. Bunches of pairs (BoPs): minimum number of 
pairs « pa i rs within maximum angular difference 

AO 

Moreover, if the pairs are due to translations in Euclidean co- 
moving space, then the vectors drawn between the two mem- 
bers of any pair should be parallel to one another. Therefore, 
a second filter can be applied, namely, selecting those sets of 
pairs where at least n pa i rs vectors connecting those pairs point 
in (nearly) the same direction. Such a collection is referred to 
as a 'bunch of pairs' (BoP). This means that any lone pair of 
objects pointing to a given direction, isolated from other pairs 
but fulfilling the redshift similarity criterion, is removed by this 
filter, so it is not taken into account at all. 

In this analysis, n pa i rs = 3 is adopted; i.e. if in a given 
direction (within some angular tolerance — see below) there 
are only two pairs of objects fulfilling the redshift similarity 
criterion and have vectors in this direction, then these pairs are 
treated as coincidences and disregarded. In practice, it turned 
out that the adoption of «p a i rs = 3 is very conservative, since 
even BoPs with «p a j rs = 5 constitute "noise". 

According to this definition, therefore, a BoP consists of 
at least 3 vectors connecting pairs of AGNs, and these vectors 
must fulfill the condition that the angle between any two of 
them is < A9, where A9 is a small angle. 

A physically reasonable value of A6 depends on possible er- 
rors in the three-dimensional spatial positions of the two AGNs 
contributing to a pair. As mentioned in Roukema (1996), an 
upper limit of a few 100 km/s on peculiar velocities limits er- 
rors in positions to about 1-10 h~ x Mpc. For pair separations 
of about 3000-6000 h~ l Mpc, this implies maximum errors of 
about 0.003 radians. However, to be conservative, in order not 
to remove too many pairs that might correspond to a true signal, 
a much larger tolerance can be adopted. 

In general, particularly given the non-uniform distribution 
of the objects with known redshifts in the catalogue we used, 
varying AO has a clear effect on the distribution of BoPs. If 
a high value of A8 is set, several BoPs are merged together, 
or "smeared", while for lower values, both the "noise" and 
a potential signal are suppressed. The smearing effect can be 
monitored by checking the distribution of the angles of vectors 
within BoPs. For the present study, a A8 which is low enough to 
avoid the "smearing" effect of BoPs was found by starting with 
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a high value and decreasing it until the "poorly collimated" 
BoPs vanished. This value, AO = 0.075 rad, is (conservatively) 
well above the limit to which angular errors are likely to occur 
due to observational effects such as peculiar velocity and the 
precision of redshift estimates. 




2 4 

redshift z 

Fig. 1. Redshift distribution of the sample used for analysis — see 
Sect-El 

2.2.3. Selection effects and how to avoid them: L se i ec 

If a topological signal is present in an AGN catalogue, then a 
PSH filtered for short lifetimes and confined to the BoPs should 
show a higher signal-to-noise spike at the length of one of the 
generators of the 3-manifold of space than in the absence of 
these filters. However, the presently available AGN catalogues 
generally include deep and narrow surveys in which many ob- 
jects are located in a small solid angle or else in a long but 
narrow angular strip of the sky. 2 It is obvious that a pair of 
(small) patches of the sky well-sampled with observations in- 
creases the probability that there will be objects of the same 
redshift in each patch and so would mimic a topological signal 
in a PSH. Effects of clustering are likely to generate higher than 
random numbers of parallel vectors, of nearly equal length, be- 
tween pairs of objects. This is a selection effect that is likely to 
lead to non-topological spikes in a PSH. 

A third filter is possible to counter this effect. While spikes 
resulting from this selection effect are caused by objects clus- 
tered in small regions of space, spikes that represent topologi- 
cal pairs should occur for objects spread throughout space. So, 
by plotting the vectors represented in a spike in combination 
with the spatial positions of the objects from which they origi- 
nate, it should be possible to judge whether the spike is caused 
by a selection effect or whether it is consistent with topologi- 
cal lensing. More quantitatively, a minimum spatial separation 
^selec between the objects generating the pairs in a bunch can 
be defined in order to exclude most of the pairs from any BoP 

2 We deliberately ignore the SDSS for this very reason — see 
Sect.[0]andEl 



generated by this selection effect. Any spike of topological ori- 
gin should remain after the application of the L se i ec criterion, 
while artefactual spikes, caused by this selection effect, should 
be removed. 




16 18 20 22 24 

B 

Fig. 2. Apparent R (upper panel) and B (lower panel) magnitude dis- 
tribution of the sample used for analysis — see Sect. 13.11 

3. Results 

3.1. The sample of AGNs 

We used an early version of the catalogue compiled by Flesch 
& Hardcastle (2004) 3 which aligns and overlays the year 2000 
releases of the ROSAT, HRI, RASS, PSPC, and WGA X-ray 
catalogues, the NVSS and FIRST radio catalogues, the Veron 
QSO catalogue, the Principal Galaxy Catalogue, and the SDSS- 
EDR onto the optical APM and USNO catalogues. The version 
we used is dated 18 October 2001. It contains 26533 objects 
with known redshifts, but we actually took 15762 with z > 0.1 
for our analysis. The redshift and apparent magnitude distribu- 
tions of this sample of 15762 objects are shown in Figs^ an d 

® 

3 Available at ftp://ftp.quasars.org 



Marecki, Roukema & Bajtlik: Cosmic crystallography using AGNs 



5 



The dominant observational uncertainties that are relevant 
for pair separation histogram analysis are the redshifts. Angular 
positions on the sky are known for these objects to better than 
an arcsecond, which corresponds to an error of much less than 
a megaparsec in the sky plane. 

Typical redshift uncertainties are ±0.01 in the red- 
shift value, corresponding to measurement errors of about 
30 hr x Mpc at z = 0.10 or about 10 \T X Mpc at z = 3.00. 
Additional uncertainty in the true distance estimates for a fixed 
set of local cosmological parameters (Q m , £2a, w, Ho) comes 
from the unknown peculiar velocities of the objects relative to 
the comoving coordinate frame. A peculiar velocity of a few 
100 km/s would give about 1-10 hr x Mpc additional radial un- 
certainty in position. 

We also attempted to explore a newer version of that 
catalogue which includes the Sloan Digital Sky Survey/Data 
Release 1 (SDSS/DR1) (Schneider et al. 2002). It turned 
out, however, that this introduced a large selection effect, as 
discussed in Sect. 12.2.31 which had to be removed anyway. 
Consequently, we abstained from trying any of the more re- 
cent but — at the time of writing — still incomplete releases of 
SDSS. 

3.2. Application of the method 

3.2.1 . Maximum redshift difference: Az 

As discussed in Sect. I2.2I a tight constraint has to be imposed 
on Az in order to remove most pairs that do not represent an 
AGN in a single active phase, while also including some toler- 
ance to reflect observational uncertainty in three-dimensional 
comoving spatial positions caused by peculiar velocities and 
inaccuracies in measured AGN redshifts. Fig.~3]shows the PSH 
for all objects we used for our analysis (Sect. I3.U assuming 
local cosmological parameters Q m = 0.3, Qa = 0.7, w = -1 
after filtering with Az/z = 0.005. 

No obvious signal is present in this figure; despite removing 
a large part of the noise, either no signal is present at all, or 
any signal present is still hidden. As the following discussion 
shows, non-topological signals are present but require further 
filtering to successively detect and then remove them. 

3.2.2. Bunches of pairs (BoPs): minimum number of 
pairs 7ip a i rs within maximum angular difference 

A9 

As described in Sect. 12.2.21 pairs caused by topological gener- 
ators should be parallel to each other. Fig.~4]shows the result of 
following the Az/z = 0.005 criterion by looking for BoPs with 
A9 = 0.075 rad and «p a j rs = 3. 

As can been seen in Fig. |ff] the concept of BoPs works 
quite well as a filter of the "noise" introduced by pairs of non- 
topological origin: several spikes with good signal-to-noise ra- 
tios clearly become visible when in addition to taking into ac- 
count AGNs' short lifetimes — an effect of the Az filter — we 
also reject all the pairs that do not constitute BoPs with at least 
three members (« pa i rs = 3). 
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100 



50 
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Fig. 3. PSH (pair separation histogram) for the objects in our sam- 
ple. The only filter applied is requiring the redshifts to both objects in 
a pair to be very similar — Sect. l2.2.T1 The error tolerance in the deter- 
mination of redshifts is 0.5% difference in redshifts in a pair. The hor- 
izontal axis shows separations of objects in pairs in units of lr [ Mpc. 
The vertical axis shows the number of pairs in a separation bin equal 
to IhT 1 Mpc. 




2000 3000 4000 5000 6000 7000 8000 



Fig. 4. PSH of Fig. [~] modified according to the BoP filter 
fSect. l2.2~2l . The horizontal axis shows separations of objects in pairs 
in h~ l Mpc. The vertical axis shows the number of nearly parallel pairs 
in a bunch. As in the case shown in Fig. [5] we allowed for 0.5% dif- 
ference of redshifts in each pair. The angular separations of pairs in a 
bunch of pairs is A8 < 0.075 rad, and « Dau - s = 3. 



While the signal-to-noise has again been improved, this is 
still insufficient to establish the topological nature of the signal. 
In particular, the apparently good signal may be the result of 
selection effects. Clearly, at least some of those spikes are likely 
to be non-topological in origin, since there is a distribution of 
many spikes close together, especially for the spikes at around 
2 - 3/T 1 Gpc, whereas the error in distances for a true spike 
should not be larger than about 10h~ l Mpc. If any of the spikes 
are topological, then they just happen to have a signal of about 
the same amplitude as noise (or selection effect) spikes. 
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Fig. 5. The distribution of directions the BoPs (Sect. 12.2.21 point 
to in galactic coordinates. The plots are in polar coordinates: the 
north galactic pole is in the centre, circles are plotted for b = 
75°, 60°, 45°, 30°, 15°, 0°. The upper panel is for pairs shorter than 
5/r 1 Gpc; the lower panel is for pairs longer than 5/T 1 Gpc. The radii 
of filled circles marking the positions are proportional to the number 
of pairs in a particular BoP. 
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Fig. 6. Locations of pairs in the BoPs represented by the two highest 
peaks around 7000/i _1 Mpc in Fig. [3] shown as their projections onto 
the equatorial plane of the cartesian coordinate system in which the 
direction of the z axis is aligned with the direction the BoP points to. 
(The coordinates are in h~ l Mpc.) An overwhelming majority of pairs 
are highly concentrated in one or two regions, which is what would 
be expected from selection effects which mimic a topological signal. 
A real topological isomorphism should generate BoPs throughout the 
plane, not just in a few regions. 



3.2.3. Selection effects and how to avoid them: L se i ec 

Before applying the L se i ec criterion to remove selection effect 
spikes (Sect. 12.2.31 . it is interesting to look at the orientation 
of all the BoPs of Fig. @] which is plotted in Fig. |5] Since 
the peaks happen to emerge in two groups (2 - 3/T 1 Gpc and 
6 - 8/T 1 Gpc), two plots are presented, one for each of these 
groups. It turns out that the bifurcation we observe in Fig. |4] 
reflects the bias caused by the Galaxy, i.e. the scarcity of ex- 
tragalactic observations near the galactic equator (the so-called 
zone of avoidance). The BoPs present in the upper panel of 



Fig. [5] where the distances in pairs are < 5hr l Gpc (the "short" 
ones), connect objects mostly on the same side of the galac- 
tic equator, either north or south. Therefore, the "short" BoPs 
predominantly point towards low galactic latitudes. 

On the other hand, the BoPs present in the lower panel 
connect the objects on both sides of the galactic equator, and 
so they largely point to high galactic latitudes. This explains 
why the histogram is bimodally distributed, but the question of 
whether any spikes are topological remains. 
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Fig. 7. PSH from Fig. |4] with the L selec = 150/T 1 Mpc criterion 
imposed. The horizontal axis shows separations of objects in pairs in 
h~ l Mpc. The vertical axis shows the number of nearly parallel pairs 
in a bunch. As in the case shown in Figs. [3] and [4] we allowed for 
0.5% difference of redshifts in each pair; and as in Fig.0 the mini- 
mum angular separation of pairs in a bunch is Ad < 0.075 rad with 
a minimum of «p au - s = 3 pairs. Clearly, the £ se ] ec criterion removes 
most selection effect BoPs. 
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Fig. 8. Locations of pairs in the BoP represented as the highest peak 
in Fig. Q shown in a cartesian coordinate system analogous to that in 
Fig. [6] Thanks to the £ se ] ec criterion, there are no selection-effect- 
generated clumps of pairs. The remaining pairs are either the tail of a 
distribution of BoPs satisfying all the filters; or they could possibly be 
a topological signal. 



The most outstanding spikes in Fig. @] are the two at about 
7000/z -1 Mpc. However, as mentioned in Sect. 12. 2.31 what we 
could and should expect is an ubiquity of false signals result- 
ing from an uneven distribution of objects with known red- 
shifts, since spectroscopic observations are often concentrated 
in small solid angles. 

A visual way to investigate this for a particular peak in the 
histogram of BoPs is to plot the distribution of pairs in a given 



BoP in the cartesian coordinate system where the z axis points 
to the averaged direction that the BoP points to. If a given BoP 
is a result of multiple connectedness, a global property of the 
Universe, then the distribution of pairs should not favour any 
particular area on the "equatorial plane" of this coordinate sys- 
tem; in other words, the isomorphism from one copy of the fun- 
damental domain to another applies equally well everywhere 
within the fundamental domain — it is a Clifford transforma- 
tion. (This is by assumption — "generator pairs", or "type II 
pairs", only give spikes in PSHs if Clifford transformations are 
present.) 

The spatial distribution of pairs in the BoPs responsible 
for the two peaks in Fig. @] around 7000/i~ 1 Mpc are shown 
in Fig. [6] Clearly, most of these nearly parallel pair vectors 
come from a selection effect in the compilation of catalogues 
we used. The plots for the other high spikes in Fig.|4]are very 
similar. 

As described in Sect. 12.2.31 this type of selection effect 
should be removable by applying the i se i ec criterion. We tried 
several values for this parameter, and found that the i se i ec cri- 
terion works effectively for values £ se i ec > 150/z 1 Mpc. By 
applying L se i ec = 150/z -1 Mpc, the "clumped" pairs in the 
BoPs represented in Fig.|4]disappear and, as a result, the L se i ec - 
filtered BoPs contain much less pairs. The histogram is now the 
one shown in Fig.0 clearly the BoPs shown in Fig.|6]no longer 
contain many more pairs than any other BoPs, of which most 
must be BoPs representing noise. 

There is no longer any obviously significant signal among 
the remaining BoPs. The highest peak represents the most 
abundant BoP, with 8 pairs. Fig. [8] shows the spatial distri- 
bution of the pairs for this particular BoP. The averaged di- 
rection they point in has the following galactic coordinates: 
/ = 241.4°, b — 78.2°. The average length of the vectors in 
this BoP (i.e. the length of the putative topology generator) is 
651 1/T 1 Mpc. 

4. Discussion, conclusions and future work 

Several filters that remove noise from the pair separation his- 
togram (PSH) have been presented and applied: these clearly 
show that the PSH method can be improved on to find weak 
signals and to remove false signals. 

The remaining best peak in the filtered PSH, from the full 
AGN compilation, is somewhat similar in its low significance 
level to that of Weatherley et al. (2003). They repeated the 
search made by Fagundes & Wichoski (1987); i.e. assuming 
that the Milky Way was once a quasar, they looked for its 
"ghost" images using up-to-date quasar catalogues including 
the 2dF Galaxy Redshift Survey (2dFGRS) 4 (Colless et al. 
2001, 2003) and SDSS/DR1 (Schneider et al. 2002). Out of 
the 7 antipodal pairs that might be potential topological images 
of the Galaxy, 6 involve objects belonging either to 2dFGRS 
or to SDSS or both. So, removal of the 2dFGRS and SDSS 
data from the analysis by Weatherley et al. (2003) would yield 
only one solution. As in our own analysis, an apparent 7-fold 



4 The 2dFGRS is now integrated with the 2dF QSO survey (Colless 
et al. 2003). 
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increase is capable of mimicking a strong topological signal; 
the authors conclude that their result shows no significant dif- 
ference in the number of pairs found over the number expected 
purely by chance. 

It seems there are two ways to correct for this inhomogene- 
ity. The first one is obvious: one has to use a comprehensive, 
large area (whole sky, at best) catalogue obtained via a uniform 
procedure. The SDSS, when completed, will be a very good re- 
source for tackling this problem, even though it will only cover 
n steradians 5 . 

Note, however, that the early SDSS releases, especially 
SDSS/DR1 (Schneider et al. 2002) 6 have an intrinsic problem 
that outweighs any of the normal observational problems such 
as apparent magnitude incompleteness. Since these observa- 
tions are carried out in relatively narrow strips (~ 5 - 10°), 
of which two of the main ones lie along opposite sides of the 
celestial equator — and one lies fairly near the celestial equator 
— a large fraction of the quasars in the catalogue necessarily 
lie in what is nearly a single plane. So, many separation vec- 
tors will necessarily be nearly parallel and have nearly identical 
lengths, because that is the geometry of the catalogue. 

As the SDSS increases to cover a large solid angle, this 
problem will become less important but will perhaps persist. 
Thus, the filters presented in this paper, A9 in Sect. 12.2.21 and 
L se i ec in Sect. 12.2.31 still will be needed. They could poten- 
tially be "tuned" to the specific geometry of the SDSS in order 
to filter out the expected "selection effect" spikes caused by 
the survey geometry. Additional filters with spatial orientations 
specifically designed for the celestial positions of the SDSS se- 
lection function, in contrast to the generic filters presented here, 
could also help remove false signals. 

It is interesting to note at this point that the SDSS has its 
radio counterpart: the VIA Faint Images of the Radio Sky at 
Twenty-cm survey (FIRST, White et al. 1997); and this FIRST 
catalogue might be used instead of SDSS. The radio oriented 
approach has an important virtue that might help not only to 
detect a possible topological signal in the future when the com- 
plete SDSS data becomes available, but also help to further 
test a candidate generator and to choose the right topological 
model, if such a signal exists. Instead of considering AGNs as 
mere point-like objects, as has always been done so far, one 
can explore their morphologies. The FIRST catalogue contains 
81 1,000 sources, which implies that at the apparent magnitude 
m(v) ~ 24 limit of SDSS, ~ 50 % of the optical counterparts to 
FIRST sources will be detected. 

The availability of the maps (two-dimensional images) of 
radio objects listed in the FIRST catalogue provides a com- 
pletely new possibility to design a statistical test, such as the 
PSH method combined with the filters presented above, by 
comparison of the morphologies of potential pairs of topologi- 
cally lensed radio-loud AGNs (RLAGNs). 

We do realise that the power of this method is limited by 
a crucial factor: radio structures of AGNs are known to evolve 
rather quickly. This requires fine accuracy in the redshift data. 
When fully available, the SDSS, with an accuracy of radial ve- 

5 See http://www.sdss.org for details. 

6 http://www.sdss.org/drl/ 



locity estimation reaching 30 km s , will give an opportunity 
to apply this method. In other words, if there is a statistically 
significant clue that a set of objects seen at identical redshifts 
could be pairs of topologically lensed images, then they should 
exhibit a substantial degree of morphological similarity, after 
correction for the projection effects caused by different view- 
ing angles. 

These projection effects imply another caveat: the beaming 
effect changes the appearance of RLAGNs considerably, so the 
search for topologically lensed images has to be confined to 
objects in nearly antipodal positions within a pair. If that con- 
straint is pushed even further, i.e. to (almost) exact antipodal 
positions, then this method becomes restricted to a special case, 
i.e. to the approach explored by Fagundes & Wichoski (1987) 
and followed-up by Weatherley et al. (2003), and provides yet 
another chance to check if we can see topological ghosts of the 
Milky Way. The difference between the two methods in that 
case is that for an RLAGN pair at a redshift z corresponding 
to distance L from the observer, Fagundes & Wichoski (1987) 
and Weatherley et al. (2003) only consider generators of length 
L, while the more general PSH method, constrained to (almost) 
exact antipodal positions, considers both generators of length L 
and of length 2L. 

Our results can be summarised as follows: 

(1) Application of the Az/z = 0.005 filter (Sect. |2~2~T1 
Sect. I3.2.T1 to a catalogue of 15762 objects with z > 0.1 
from an observational compilation of AGN catalogues en- 
ables removal of most of the roughly 10 8 pairs from the pair 
separation histogram (PSH), which for physical reasons — 
the short lifetimes of AGNs — are much more likely to be 
noise than pairs satisfying the filter. This is shown in Fig. [5] 

(2) Application of the A9 = 0.075 rad and «p a j rs = 3 fil- 
ters, which select out only those bunches of pairs (BoPs) 
which are consistent with being caused by an isomorphism 
(generator) between copies of the fundamental polyhedron 
(Sect. l2.2~2l Sect. 13.2.21 . and thereby further removes noise 
from the PSH, results in further improvement in the signal- 
to-noise, and yields many BoPs which become obviously 
visible — see Fig.|4] 

(3) However, these BoPs are best interpreted as selection effect 
spikes (Sect. lzZOl . A length criterion L se j ec 150 hr x Mpc 
to remove those BoPs most likely to be selection effects is 
effective (Sect. l3~2"""t in doing so. 

(4) The highest spike remaining after these filters (Fig. [7), 
points to the galactic direction (/ = 241.4°, b - 78.2°) and 
is of length 65 1 \hr l Mpc. However, it is not highly signifi- 
cant. 

When the full SDSS becomes available, applying this 
method would definitely be worth trying in order to extract 
any signal that might be present. Even though WMAP anal- 
yses have so far indicated that no topological signal is likely to 
be present at scales much below the horizon scale, it would be 
imprudent to rely on one observational method alone to con- 
sider such a fundamental cosmological property — the shape 
of the Universe — to be unmeasurable on these scales. 

Another prospect for future work is that if an interesting 
signal were found, then comparing radio maps of the putative 
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topologically lensed RLAGNs would seem to be a promising 
tool for discriminating between true and false topological sig- 
nals. 
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